TUM/T39-98-6 

On the evaluation of some three-body variational integrals 



Jose Caro 

Physik Department, Technische Universitat-Miinchen, D-8 5747- G arching, Germany 

Abstract 

Stable recurrence relations are presented for the numerical computation of 
the Calais-Lowdin integrals 

J dr 1 dr 2 r l 1 ~ 1 r^ 1 r^ 2 1 exp {-an - fir 2 - jr 12 } 

(I, m and n integer, a, (5 and 7 real) when the indices I, m or n are negative. 
Useful formulas are given for particular values of the parameters a, (3 and 7. 
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I. INTRODUCTION 



When dealing with the three-body variational problem with Hylleraas basis, it is usually 
necessary to make extensive use of integrals of the general form |l] 

/(/, m, n; a, ft, 7) = J dr 1 dr 2 r^V™ -1 ^" 1 exp {-an - [3r 2 - 7^12} , (1) 

where n = |ri|, r 2 = 1 1*2 1 and ri% — \?2 — ri|. 

For the case of /, m and n non-negative (that is, non-negative powers of ri, r 2 and 
ri2 once the volume element has been taken into account), powerful, simple and stable 
recurrence relations that permit the numerical calculation of these integrals can be found in 
the literature 0. However, it is sometimes essential to have also an expression for one of 
the integer indices being negative. For instance, that happens in the atomic problem when 
one wants to consider the mean value of the rf 2 2 operator [[3] or relativistic corrections [Q ; or 
in the nuclear problem when non-local terms are included in a Yukawa-like interaction ||. 
In some cases, the integrals must be computed in every step of the non-linear optimization 
procedure, and hence it is clear the need of having a quick and reliable algorithm to compute 
them. The specific cases 1(1,1,-1) and 1(0,-1,-1) were already considered in Refs. || 
and || respectively. For 7 = much work has been done PJTHTTH, also including explicitly 
the coupling of the angular momentum of the two dynamical particles ||12|| . Some work has 
been devoted to the analogous integrals for four- or more-body problems ||0.[TT].[l^.[ri . 

The method proposed in this work to obtain the integrals (|l|) is specially useful when 
the same exponential coefficients a, [3 and 7 appear in several elements of the variational 
basis. 



II. GENERAL PROPERTIES OF 1(1, m, n) 

To study the general properties of the integral ([1]) for /, m and n (possibly negative) 
integer numbers and a, (3 and 7 real it is convenient to make use of perimetric coordinates 



u = —r\ + r 2 + ri2 

v = T\ - r 2 + r 12 \ (2) 
w — r 1 + r 2 - r 12 

in terms of which the initial integral reads 

1(1, m, n; a, f3, 7) = 2~^^ Ip ^ m> „. P+l, £+2, 2+£\ , (3) 

where 

/*oo /*oo /*oo 

I p (l, m, n; a, b, c) = / du dv dw (v + w) l (u + w) m (u + v) n exp {— au — bv — cw} . 
Jo Jo Jo 

(4) 
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The integral I p is explicitly invariant under permutation of conjugated pairs of parameters 
{(I, a), (m, b), (n, c)}, and therefore 

1(1, m, n; a, (3, 7) = I(m, I, n; [3, a, 7) = I(n, m, I; 7, ft, a) , (5) 

symmetry that will be used throughout this work. 

The long range convergence of I p is ensured if a, b and c are positive real numbers, that 
is, if 

a + (3 > , a + 7>0and/3 + 7>0. (6) 

That means that one of the exponentials parameters, a, (3 or 7, can be zero or negative, 
provided that the other two are bigger than the absolute value of the former. Note also 
that one of the exponential coefficients of I p can be zero if the power of the corresponding 
integration variable is negative and high enough. For instance, a = with I = and 
m = n = — 1 would yield a convergent result. Anyhow, this is an almost useless case for 
the variational problem, because for higher power integrals (that very likely should also be 
considered) a = would lead to divergent quantities. From now on, we assume that the 
requirements @ are fulfilled. 

The study of the short range convergence can be straightforwardly done case by case. 
Summarizing, for I, m and n integer, and a, (3 and 7 real such that the conditions (H) are 
fulfilled, the integral (|I|) is convergent if and only if 

1>-1, m> -I. n > -1 and I + m + n > -2 . (7) 

To have a procedure to generate the whole set of integrals (HD one needs relations for the 
cases 1(1, m, —1) and /(/, —1, —1) where I and m are non-negative. 

As soon as we have checked that the integral we are looking for is convergent, integration 
over one parameter can be applied to lower the conjugated power, 

/■oo 

1(1, m,n; a, /3, 7) = / dcl(l,m,n+ I; a, (3,c) . (8) 
h 

On the other hand, derivation can always be used to increase indices, 

(-d a ) p I(l, m, n; a, f3, 7) = 1(1 + p, m, n; a, f3, 7) . (9) 

These properties, together with 

1(0, 0, 0; a, (3, 7) = (a + (3)- l (a + ^((3 + 7)- 1 , (10) 

are useful to derive all the integrals. Note also that for A > 

/(/, m, n; Xa, \(3, A7) = A^ +m+n+3 )/(/, m, n; a, (3, 7) , (11) 

that is, for given /, m and n, I is a homogeneous function of a, (3 and 7. This fact, together 
with properties (|j) and @ yields a quite general recursion. Indeed, differentiating with 
respect to A in the equation above one gets the recurrence relation 
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(I + m + n + 3)1(1, m, n; a, /3, 7) = al(l + 1, m, n; a, (3, 7) + (31(1, m + 1, n; a, (3, 7) 

+ jl(l,m,n + l;a,0,j), (12) 



valid for well defined integrals, in our case /, m and n fulfilling conditions (|7]). In general, 
this recursion is of little utility, for to use it downwards, which is the obvious direction, one 
would have to know the value of the integrals on a plane l + m + n = constant. We will take 
profit of a particular case of Eq. ( |i~2"|) in section |TV|. 



III. CASE 1(1, m, -1) WITH l,m > 

For the family of integrals 1(1, m, — 1; a, ft, 7) with l,m > 0, a variation of the method 
exposed in Ref. |§ can be applied. The recurrence relation that one gets is the following, 

1(1, m, -1; a, f3, 7) = — !— [11(1 - 1, m, -1) + ml(l, m - 1, -1) + B(l, m)] , (13) 

a + p 

where 

/•OO 

B(l,m;a,$,i) = l\m\ / dc (a + c)"'" 1 (/3 + c)"™" 1 , (14) 

which is a symmetric function under (/, a) ^ (m, (3) exchange, can be obtained through the 
relation 

B(l, m) = — — [IB (I - 1, m) - mB(l, m - 1) + C(l, m)\ . (15) 
a — [3 

Here the function C(l,m) reads 



C(l,m;a,(3,i) 



' (m - l)\ ((3 + ~f)- m if / = and m>0 

-(/ - 1)! (a + 7)"' if I>0 and m = 

log{g^ if £ = and m = (16) 

otherwise 



and is defined so that the recursion flTo] ) holds also for / = m = although -8(0, — 1) and 
B(— 1,0) are divergent. Note that C(l,m) is antisymmetric under (/,a) <->• (m,f3). 

Unfortunately, in the recursion (|15[) subtractions are involved, and hence one must look 
over the stability against roundoff, in particular when a and (3 are close to each other. 

It is also possible to relate B(l,m) to Gauss hypergeometric function, 2 Fi [0, yielding 

/! 777 1 

B(l, m; a, [3, 7) = —— (a + 7)^ (/? + iT™ 2^(1, / + 1; m + / + 2; z) , (17) 
m + £ + 1 

where z = ■ The use of the integral representation of the hypergeometric function gives 

r i f l (1 — f) m 

B(l, m; a, (3, 7) = (I + m)\ (a + ((3 + 7 )" m J dt ± _ J . (18) 
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From the definition fll4|) it is possible to prove the equation 

B(l + l,m) + B{l,m + l) = l\m\{a + -f)- {l+1) ([3 + ~f)- (rn+1) (Z,m>0). (19) 
Plugging this relation in Eq. (|i"5|) yields 



(Z + m) B(l - 1, m) - (a - /3) B(Z, m) - (Z - 1)! m! (a + 7 )~' + 7)" 







(20) 



valid for m > and Z > 0. This equation permits to lower one unit the index I of B(l,m) 
with numerical stability if a > (3. In the opposite case, the symmetry of B(l,m) can be 
used to lower the index m (see Fig. p. 

On the other hand, using the Gauss relations for contiguous hypergeometric functions 
one obtains 



mB(l + l,m-l) + (m -l£)B(l,m) - l£B(l - l,m+ 1) = , 



(21) 



where £ = 1 — z = ^p-. This relation defines a recursion that can be used to move on the 
diagonals m + I = constant. As it is shown in Fig. [I], the straight line m(Z) = £Z on the 
Z-m-plane separates the stability regions of the recursion (PH), so that one can move with 
stability from this line in diagonal steps. 

The final recipe to compute the set of 1(1, m, — 1; a, f3, 7) for l,m < N is the following 
(see Fig. [I]). First, two 5's are to be computed numerically to the required accuracy, namely 



B 



2N 
14* 



,2iV 



27V 
14£ 



and -B f 



2N 
14£ 



+ 1,27V 



2V 



1+? 



- 1 



(22) 



(respectively, points Pi and P 2 in Fig. [T]). Then the recursion fl2T|) is used to generate all 
needed starting points to use the recursion ( p0|) leftwards (downwards) if a > (3 (a < (3). 
Finally, the B J s obtained in this way are introduced in Eq. ([13]). To generate the two 
initial B's one can compute the integral in Eq. flT8| ) by Gauss- Legendre quadrature. To 
optimize the computation of the quadrature a change of variable is needed. First, we use 
the symmetry of B(l,m) to render < z < 1. Next, we apply in Eq. fll8|) the change of 
variable (t — > s = s(t)) 



S ^ { log (2 - z - t) 



if < z < 0.8 
if 0.8 < z < 0.99 



(23) 



For values of z greater than 0.99 the hypergeometric function can be computed using the 
transformation formula 15.3.11 of the Ref. fl6| . With the prescription above, more than 



fifteen stable figures were obtained using 32 Gauss-Legendre points for 2N < 60. Note that 
the prescription ( p3| ) has been optimized for the computation of the two initial 5's given in 
the expression (p2|) , and will not provide a similar accuracy for arbitrary values of I and m. 

The particular case a = j3 is specially simple. Indeed, in that case the B function to be 
included in Eq. (|13|) is 



B(l, m; a, a, 7) 



l\m\ 



I + m 



.( a + 7 )-G4m+l) 



(24) 



and the calculations are numerically stable. The case a = (3 is not only a mere academic 
example. In many practical problems the variational basis is chosen so that any element has 
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the same exponential coefficient both for the coordinates r\ and r^- If the physical problem 
requires to deal with 1(1, m, —1) integrals, then it is sensible to check whether such a basis 
can produce the required accuracy. This selection was successfully used in the context of a 
nuclear theory problem ||. 

In Table [I] we give some particular values of /(/, m, —1; a, (3, 7) with fourteen significant 
figures to provide the reader with checking points. 



IV. CASE 1(1, -1, -1) WITH I > 



To generate the set of integrals /(/, —1, —1) use can be made of the relation 

a 1(1 + 1,-1, -1; a, (3, 7) = (/ + 1)1(1, -1, -1; a, (3, 7) 

- (31(1, 0, -1; a, (3, 7) - -yl(l, -1, 0; a, (3, 7) , 



(25) 



which is valid for I > 0. This equation is easily obtained as a particular case of recurrence 
(p"2|). For I = direct calculation yields 



I(0,-l,-l;a,(3 n ) = — 

la 



vr 2 a + 7, a + (3 T . / (3 - a 
T ~ lo S IT, — lo S IT, Ll 2 IT, — 



- Li 2 



7 — a 



(26) 



(see also Ref. ||), where Li 2 (z) = — Jq dy y _1 log(l — y) is the dilogarithm function. 

As said, the expression (^) is not applicable for the case I — m — n — — 1. Instead, one 
gets 



7T 



al(0, -1, -1; a, (3, 7) + (3I(-1, 0, -1; a, (3, 7) + 7 /(-l, -1, 0; a, (3, 7) = — . (27) 
To fix the constant in the right hand side, we had to make explicit use of the expression 



The recursion (25), which in general is numerically unstable upwards, can be used with 
stability to decrease the index / if a > 0, which is the interesting case in physics. But then, 
one needs as starting point the integral with the highest wanted /. As it can be derived from 
Eqs. (||-(T0]), that integral can be obtained through the computation of the quadrature 



/G,-i,-i;«,A7) = § (-^ 



l+l 



r 1 1 

Jo t 



2a 



2(3 



/3 + 7'/3 + 7 ' t 



(28) 



where we have defined 
1 



Gi(a,b; y) 



1 



+ E 

m=l 



m 



a + y) 

g + y 

a + b 



—1 lo § 



(a + b + 2(y-l))(a-b + 2y) 



+ 



a + y 



(a + b)(a-b + 2) 



a + y 



a-b + 2 



a + y 



b + 2y y 



a + b + 2(y-l) 



(29) 



Note that the integrand is positive, and that the sum in the function Gi is very efficiently 
computed in a single loop. For values of a, (3 and 7 of the same order of magnitude the 
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quadrature converges very quickly for not very small values of I (I > 5). This is not the case 
when one of the parameters is larger than the other, but then a simple change of variable 
helps to recover convergence. For instance, the following prescription of changes of variable 
(i-s = s(t)) 

S (f) = l M t+ f\ , if a> iW + T) (30) 
1 lo S V 1 + a+ 2 7/3+7)' 7} ~ *) otherwise 

produces the accuracies shown in Table |I|. For I = 10 (respectively, 20, 30 and 40) more than 
3800 integrals per second (respectively, 2500, 1800 and 1400) were obtained with more than 
fourteen stable figures (32 Gauss- Legendre points, see Table |TT]) in an inexpensive computer 
(a PC with a 200 MHz processor). 

The prescription ( |30| ) is not useful if both /3 and 7 are much smaller than a (e.g., 
(3, 7 < 0.01a). However, the integrals for this case can be safely generated without significant 
loss of accuracy using the recursion ( p5|) upwards. The relative error in the l-th integral 
accumulated because of cancellations, which grows with I, can be then approximated as 

„ /w , 1 . 7r 2 l\ machine-precision . , 

1(7(1,-1,-1)) = -^ . (31) 

For example, 1(60,-1,-1; 1,0.01,0.01) __ 0.768 ■ 60!, and the relative error due to cancel- 
lations in the repeatedly use of fl2"5] ) to obtain this integral is only of about three times the 
machine-precision. For smaller values of j3 and 7 the accuracy is bigger. Note that Eq. 
( p6f) is not appropriate to evaluate 7(0, —1, —1) when j3 and 7 are smaller than a. A better 
expression for this case is 



7(0,-1,-1; a, /3, 7 ) = -^ 



\a + (3J \a + -f / 



. a + (5 . (a + (5) (a + 7) , a + 7 , (a + (3) (a + 7) 
log log log log 



a '" (a — 7) 2 a ' (a — (3) 2 

(32) 



+ 2 1og^±2log (a + ' 3)(tt + 7) 



a (a — /3) (a — 7) 

where the dilogarithm function can be computed through the expansion L^x) = Z)fc_=i fr- 
For x < 0.02, corresponding to (3, 7 < O.OIck, eight terms in this expansion are enough to 
obtain sixteen stable figures. 

A few particular cases of 1(1, —1, —1) are readily obtained from the recursion (p5|). In- 
deed, for a = and / > we have 

7(Z, -1, -1; 0, P, 7) = j^j [(31(1, 0, -1; 0, (3, 7) + 7 I(Z, -1,0; 0, (3, 7)] • (33) 
The specific case a = (3 = 7 reads 

1(1, -1,-1; a, a, a) = ^ 2a y+i Sl ' ( 34 ) 



7 



where the coefficients 



S, 



oo /l+m 

EE 

m=0 \i=0 



1 



i + ll l + m + 1 



(35) 



are to be computed only once. For I < 100 one does not need more than 52 terms to achieve 
sixteen stable figures in Si without using any numerical procedure to accelerate convergence. 
The first of these coefficients is So — \- Finally, for /(/, —1, — 1; a, 0, 0) one has 



7T 2 V 

/(Z, -1,-1; a, 0,0) = —. 



4 a l+1 ' 



(36) 



Some values of 1(1, —1, — 1; a, (3, 7) are presented in Table |Tl. 



V. SUMMARY 

Some recurrence relations to compute the integrals (|l|) for all negative integer param- 
eters (I, m and n) have been presented. The stability of these recursions has been in- 
vestigated, and algorithms have been given to use them without loss of accuracy due to 
cancellations. The integrals 1(1, m, —I; a, (3, (where l,m > 0) can be generated at low 
computing cost. For the integrals /(/, —1, — 1; a, (3, 7) a quadrature involving + 1 terms is 
needed, where iV is the highest required I and a is assumed to be positive. Specially simple 
algorithms are given for the cases 1(1, m, — 1; a, a, 7), 1(1, —1, — 1; 0, (3, 7), /(/, —1, — 1; a, a, a) 
and 1(1, —1,-1; a, 0, 0). 
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TABLES 



(At) 




I = 10, m = 







1(1, m, - 
I = 10, m = 


1; 1.0,7) 
5 




1 = 20, m = 15 






I = 40, m = 


10 


(0.05,0.05) 


2 


5097803893512 


(+07) 


i 


1550294761619 


(+14) 


5 


8496691184166 (4 


-45) 


3 


7263930880406 


(+65) 


(0.05,0.20) 


1 


0069028364788 


(+07) 


l 


9856655589152 


(+12) 


1 


1263701557234 (4 


-41) 


1 


5245149068463 


(+64) 


(0.05,1.00) 


2 


1231968038912 


(+06) 


6 


1362890397889 


(+11) 


1 


8814741050023 (4 


-39) 


2 


5445500764172 


(+63) 


(0.05,2.00) 


1 


0610620020573 


(+06) 


3 


0134288353911 


(+11) 


9 


0355811977946 (4 


-38) 


1 


2652815721954 


(+63) 


(0.05,5.00) 


4 


2435098606331 


(+05) 


1 


1992798691243 


(+11) 


3 


5754444925655 (4 


-38) 


5 


0533341605295 


(+62) 


(0.20,0.05) 


5 


9511548712149 


(+06) 


1 


4149826065364 


(+12) 


7 


4333291536966 (4 


-39) 


4 


6331742767202 


(+61) 


(0.20,0.20) 


2 


5245340864963 


(+06) 


3 


8096117209107 


(+11) 


1 


5683171148516 (4 


-38) 


1 


3617534581186 


(+61) 


(0.20,1.00) 


4 


9479973346792 


(+05) 


7 


1119098709184 


(+10) 


1 


5019778106939 (4 


-37) 


2 


7821828122703 


(+60) 


(0.20,2.00) 


2 


4501414762596 


(+05) 


3 


5414728355891 


(+10) 


7 


3432360158513 (4 


-36) 


1 


3921604053583 


(+60) 


(0.20,5.00) 


9 


7730721646852 


(+04) 


1 


4149098737738 


(+10) 


2 


9191066335088 (4 


-36) 


5 


5698482446321 


(+59) 


(1.00,0.05) 


2 


3506936424270 


(+05) 


2 


1738100662595 


(+08) 


1 


8638785541149 (4 


-30) 


1 


3896745377911 


(+53) 


(1.00,0.20) 


6 


8960382056483 


(+04) 


9 


7511110023298 


(+07) 


7 


0043082969132 (4 


-29) 


2 


2092326515384 


(+51) 


(1.00,1.00) 


2 


6754225852273 


(+03) 


2 


0193624551498 


(+07) 


1 


4992836690734 (4 


-29) 


1 


5383073655665 


(+49) 


(1.00,2.00) 


9 


7180274971942 


(+02) 


1 


0013416984929 


(+07) 


7 


5130229486046 (4 


-28) 


6 


9573280539039 


(+48) 


(1.00,5.00) 


3 


5921654332679 


(+02) 


3 


9932069716804 


(+06) 


3 


0070041814429 (4 


-28) 


2 


7138972320215 


(+48) 


(2.00,0.05) 


5 


6374109387576 


(+04) 


1 


4250979709052 


(+06) 


9 


1546004403515 (4 


-24) 


2 


9503794635963 


(+49) 


(2.00,0.20) 


1 


5196492752737 


(+04) 


1 


7607751878871 


(+05) 


1 


4095247283656 (4 


-24) 


1 


9130948733180 


(+47) 


(2.00,1.00) 


1 


5320840275499 


(+02) 


3 


8998899902794 


(+04) 


7 


9782256684551 (4 


-22) 


8 


8675800592916 


(+40) 


(2.00,2.00) 


1 


6829755898961 


(+01) 


1 


6286558844196 


(+04) 


3 


5727727922307 (4 


-22) 


8 


6310437705915 


(+39) 


(2.00,5.00) 


4 


3389324770986 


(+00) 


6 


1452543996272 


(+03) 


1 


3850914702513 (4 


-22) 


2 


9003763096969 


(+39) 


(5.00,0.05) 


8 


9272368511676 


(+03) 


2 


0896049079214 


(+03) 


1 


2958485756621 (4 


-18) 


1 


7703597474580 


(+44) 


(5.00,0.20) 


2 


3569834307325 


(+03) 


5 


6914967826940 


(+02) 


1 


0121437234633 (4 


-17) 


2 


4003665222216 


(+42) 


(5.00,1.00) 


1 


5138128168377 


(+01) 


6 


1945864757004 


(+00) 


2 


8220667212483 (4 


-13) 


8 


7853148834688 


(+33) 


(5.00,2.00) 


3 


2220328921320 


(-01) 


5 


5803922147916 


(-01) 


1 


0214473659965 (4 


-12) 


1 


1612978333117 


(+28) 


(5.00,5.00) 


3 


4970794973642 


(-03) 


1 


0465458754289 


(-01) 


2 


1723867163698 (4 


-11) 


1 


6228394644288 


(+24) 


TABLE 


I. Some 


values 


for the integrals 


1(1, m, -l;l,/3,7). 




For instance, 



1(20, 15, -1; 1, 0.2, 5) = 2.9191066335088 • 10 36 . 



1 


7V GI =10 


A r GL=16 


A r GL=24 


N GL = 


5 


5.9 


8.9 


11.1 


12.6 


6 


6.1 


9,1 


12.1 


13.9 


7 


5.6 


9.6 


12.9 


14.9 


8 


5.6 


9.7 


13.6 


15.4 


9 


5.9 


10.3 


14.1 


15.3 


10 


5.4 


9.3 


14.4 


15.3 


15 


5.1 


8.7 


14.0 


15.2 


20 


5.2 


8.3 


13.6 


15.0 


25 


4.7 


8.2 


12.9 


14.9 


30 


5.0 


8.1 


12.1 


14.8 


35 


4.7 


7.5 


11.4 


14.8 


40 


4.5 


7.0 


11.0 


14.6 



TABLE II. Number of stable figures (— log 10 (relative-error)) in 1(1, — 1, — 1; a, /3, 7) obtained 
using different numbers of Gauss-Legendre points (Aql) in the quadrature (|28|) with the prescrip- 
tion @. The ratio | moved in the range [0.01, 10 4 ], and g in [10~ 4 , 10 4 ]. 



10 



J(Z, -1,-1; 1,0,7) 



( o 


01 





01) 


6 


6191662489867 


(+06) 


3 


6251447231256 


(+18) 


3 


2983926561168 


(+32) 


8 


5733728363337 


(+47) 


( o 


01 





20) 


1 


7540859652228 


(+06) 


5 


3112685642066 


(+17) 


3 


3942282862188 


(+31) 


6 


9857465097091 


(+46) 


( o 


01 





50) 


6 


8609259644037 


(+05) 


2 


0155302235322 


(+17) 


1 


3180581225243 


(+31) 


2 


7475204208958 


(+46) 


( o 


01 


1 


00) 


3 


3243271983196 


(+05) 


9 


9937617391729 


(+16) 


6 


5672190824270 


(+30) 


1 


3709460938605 


(+46) 


( o 


01 


2 


00) 


1 


6469713752668 


(+05) 


4 


9876597990304 


(+16) 


3 


2808597542740 


(+30) 


6 


8512909302258 


(+45) 


( o 


01 


5 


00) 


6 


5729176667975 


(+04) 


1 


9940611043105 


(+16) 


1 


3120390818053 


(+30) 


2 


7401331921302 


(+45) 


( o 


01 


10 


00) 


3 


2854418466598 


(+04) 


9 


9695938948872 


(+15) 


6 


5599782224617 


(+29) 


1 


3700392619726 


(+45) 


( o 


20 





20) 


4 


5604766596349 


(+05) 


3 


5127689695513 


(+16) 


4 


9063595251091 


(+29) 


2 


0412224928657 


(+44) 


( o 


20 





50) 


1 


4318580613654 


(+05) 


7 


2132201828707 


(+15) 


8 


1925006515660 


(+28) 


3 


0132668197157 


(+43) 


( o 


20 


1 


00) 


6 


1814890074087 


(+04) 


3 


2566277622792 


(+15) 


3 


8027988040087 


(+28) 


1 


4138666122607 


(+43) 


( o 


20 


2 


00) 


2 


9651034471161 


(+04) 


1 


5963963424310 


(+15) 


1 


8718211282572 


(+28) 


6 


9707881275722 


(+42) 


( o 


20 


5 


00) 


1 


1742458080974 


(+04) 


6 


3522447488220 


(+14) 


7 


4555535110623 


(+27) 


2 


7776640095065 


(+42) 


( o 


20 


10 


00) 


5 


8633250910678 


(+03) 


3 


1737770991220 


(+14) 


3 


7255334615001 


(+27) 


1 


3880777482000 


(+42) 


( o 


50 





50) 


2 


7602780782403 


(+04) 


2 


1429404549106 


(+14) 


3 


1095312824899 


(+26) 


1 


3602027847527 


(+40) 


( o 


50 


1 


00) 


7 


8596287095799 


(+03) 


4 


2261129342848 


(+13) 


5 


2077776493094 


(+25) 


2 


0657542275279 


(+39) 


( o 


50 


2 


00) 


3 


2906652367190 


(+03) 


1 


8845025819044 


(+13) 


2 


3675994785514 


(+25) 


9 


4583791413986 


(+38) 


( o 


50 


5 


00) 


1 


2667740672550 


(+03) 


7 


3501151392020 


(+12) 


9 


2602020503771 


(+24) 


3 


7039195719454 


(+38) 


( o 


50 


10 


00) 


6 


3029915005903 


(+02) 


3 


6623979936588 


(+12) 


4 


6157947923024 


(+24) 


1 


8465337659004 


(+38) 


( 1 


00 


1 


00) 


9 


2233071085062 


(+02) 


3 


9012855025163 


(+11) 


3 


1367380768282 


(+22) 


7 


6500307975383 


(+34) 


( 1 


00 


2 


00) 


2 


1119403628357 


(+02) 


6 


5569071489222 


(+10) 


4 


6001747109760 


(+21) 


1 


0314129070250 


(+34) 


( 1 


00 


5 


00) 


7 


2422587486309 


(+01) 


2 


3596051913831 


(+10) 


1 


6728442527147 


(+21) 


3 


7666888585022 


(+33) 


( 1 


00 


10 


00) 


3 


5623476884278 


(+01) 


1 


1649070967644 


(+10) 


8 


2662607382138 


(+20) 


1 


8620745501178 


(+33) 


( 2 


00 


2 


00) 


8 


9219060942223 


(+00) 


6 


4071580456640 


(+07) 


8 


8443126974508 


(+16) 


3 


7177770572299 


(+27) 


( 2 


00 


5 


00) 


1 


3397156572731 


(+00) 


7 


4723481641427 


(+06) 


9 


1598047862917 


(+15) 


3 


5719240038509 


(+26) 


( 2 


00 


10 


00) 


6 


2645017020798 


(-01) 


3 


5448536509329 


(+06) 


4 


3596533454989 


(+15) 


1 


7025891413660 


(+26) 


( 5 


00 


5 


00) 


3 


7697615446040 


(-03) 


2 


6042296175688 


(+01) 


3 


4861019106182 


(+07) 


1 


4250284687126 


(+15) 


( 5 


00 


10 


00) 


6 


8653568463440 


(-04) 


3 


7115610730907 


(+00) 


4 


4359819939145 


(+06) 


1 


6880938843748 


(+14) 


(10 


00 


10 


00) 


4 


5242997355095 


(-06) 


7 


2459128782254 


(-05) 


2 


2554462512897 


(-01) 


2 


1460242640172 


(+04) 



TABLE III. Some values for the integrals 
7(10, -1, -1; 1, 0.01, 10) = 3.2854418466598 • 10 4 . 



/(*, -1,-1; 1,0,7)- 



For instance, 
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FIGURES 



(N,N) 



absolute stability 
stability for a>(3 
starting point 




FIG. 1. Stability lines of the recursions for the calculation of 1(1, m,— 1). The solid arrows 
refer to the stable flux of the recursion (21). The open ones refer to the recursion (|20|), but only 
if a > f3. In the opposite case, the symmetric of Eq. (20) under (I, a) <-> (m,(3) exchange can be 
used to move downwards with stability. 
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